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Abstract 

An important parameter to characterize the scattering matrix S for quantum- 
chaotic scattering is the width r corr of the S-matrix autocorrelation function. 
We show that the "Weisskopf estimate" dj (2ir) Tc (where d is the mean res- 
onance spacing, T c with < T c < 1 the "transmission coefficient" in channel 
c and where the sum runs over all channels) provides a good approximation to 
Tcorr even when the number of channels is small. That same conclusion applies 
also to the cross-section correlation function. 



1. Purpose 

Quantum-chaotic scattering is an ubiquitous phenomenon. It emerges when- 
ever Schrodinger waves are scattered by a system with chaotic intrinsic dynam- 
ics. Examples are the passage of electrons through disordered mesoscopic sam- 
ples, and compound-nucleus scattering. Moreover, it occurs when electromag- 
netic waves of sufficiently low frequency are transmitted through a microwave 
cavity with the shape of a classically chaotic billiard. In all these cases, chaotic 
scattering is due to the numerous quasibound states of the system that appear 
as resonances in the scattering process and that obey random-matrix statistics. 

The generic approach to quantum-chaotic scattering [l[ is based upon a 
random-matrix model for the resonances and, thus, for the scattering matrix 
S a b{E), a function of energy E where a, b denote the open channels. Within that 
approach, the energy correlation function of the scattering matrix (the ensemble 
average (S ab (E — e/2)S* d {E + e/2))) can be worked out analytically Q as a 
function of the energy difference e, and approximate expressions for the cross- 
section correlation function are also available 0, 0, 0] ■ The correlation width of 
the cross section turns out to be rather close to that of the scattering matrix 
in all cases That is why we focus attention on the S'-matrix correlation 
function in what follows. 

The analytical expression for the S'-matrix correlation function is given in 
terms of a threefold integral Q. The numerical evaluation of that integral is 
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quite cumbersome. To gain an orientation of what to expect in a given situation, 
a simple approximate expression for the width r corr of the S-matrix correlation 
function (and, by implication, of the cross section) would therefore be helpful. 
A commonly used approximation applicable in the regime of strongly overlap- 
ping resonances is the "Weisskopf estimate" . It has for example been applied 
to resonance spectra obtained from microwave experiments on quantum chaotic 
scattering (H], 0] ■ The measurements were performed in the regimes of isolated 
and weakly overlapping resonances and the associated S matrix comprised two 
dominant scattering channels and a large number of weakly coupled ones. A 
motivation for the present paper was to test the accuracy of the Weisskopf es- 
timate under such conditions. We will demonstrate that it provides a good 
approximation for r corr not only in the regime of strongly overlapping reso- 
nances. For simplicity and without loss of generality we confine ourselves to the 
case where the average S'-matrix is diagonal, (S a b) — (S a a)S a b- The unitarity 
deficit of the average S'-matrix is then measured by the transmission coefficients 
T c = 1 - |(S CC )| 2 . These obey < T c < 1 for all c. 

Naively, one might consider two alternatives for estimating r corr . (i) The 
Weisskopf estimate expresses the total average resonance width in terms of the 
mean resonance spacing d and of the transmission coefficients T c , 

c 

The sum in Eq. (flj runs over the open channels. 

(ii) The "Moldauer-Simonius sum rule" gives the following expression 

for the mean distance (1/2) (T^) of the poles of the scattering matrix (labeled 
by a running index fj,) from the real energy axis. 

( r M) = -^EHi-r c ]. (2) 

c 

For the case of unitary symmetry, the sum rule Eq. © has been derived rigor- 
ously There is no reason to doubt that the sum rule Eq. ([2]) holds also in 
the orthogonal case although a proof exists only in fragmentary form [lH . 

The width Tw in Eq. (fTJ and the double average pole distance (r M ) as given 
by Eq. ([2]) agree whenever T c <C 1 for all c. In general, however, the values of 
both quantities differ widely. For instance, for the case of a single channel with 
T « 1, Eq. (JD yields T w « d/(2%) while Eq. © yields (T^) » (1/(2%). An 
identification of Tw (of (r M )) with the correlation width r corr would suggest 
that we deal with isolated (with strongly overlapping) resonances, respectively. 
It is known 12 , 13|, E3] that Eq. ([2]) fails when any of the T c is close to unity, 



and a comparison of the values of F corr given in the figures below with Eq. © 
confirms that fact. We ascribe that failure of the Moldauer-Simonius sum rule 
to the fact that the fluctuation properties of the scattering matrix depend not 
only on the location of the poles of S but also on the values of the residues. 
Little is actually known about the latter 15, 16j |. 
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That leaves us with Eq. ([T]) as the only viable alternative. We recall the 
conditions under which Eq. (pj is obtained [f|. One uses a time-dependent 
description and considers a scattering system with constant resonance spacing 
do coupled to a number of channels. The frequency with which a typical wave 
function of the system approaches the entrance of a given channel c is do/h, 
the probability with which the system escapes into that channel is given by T c , 
the partial width for decay into channel c is accordingly d$T c / (2tt). Summing 
over all channels and postulating that the result applies also to systems that do 
not have a constant resonance spacing do, one replaces do by the actual mean 
resonance spacing d and arrives at Eq. (JXJ) . The argument being semiclassical, 
one expects Eq. (JXJ) to give an approximate expression for the average resonance 
decay width in the case of many channels or, more precisely, for JD C T C ^ 1. 

That argument leaves open the question how T\v relates to the correlation 
width r corr . In Ericson's work [17[ the identity of T\v and of r corr was postulated 
for XL-T C A proof for that assertion became available with the work of 



Ref. [18j. There it was shown that an expansion of the ^-matrix correlation 



function derived in Ref. Q in inverse powers of J] c T c yields as the leading term 
a Lorentzian with width This result implies r corr = Tw in the Ericson 

regime of strongly overlapping resonances, i.e., for J^ C T C ^> 1. A second case is 
that of many open channels each coupled weakly to the resonances. Using the 
analytical expression for the S'-matrix correlation function Q Harney et al. (l9| 
showed that in that case, T\v also provides a good approximation for r corr . 
Apart from these results for the regimes of strongly overlapping and of isolated 
resonances coupled to many channels no simple analytical expression exists for 
r colT . In the present paper we investigate how much T-yv and r corr differ for 
general values of the number of channels and of the coupling strength in each 
channel. We use the analytical result of Ref. Q for the S'-matrix correlation 
function to compute the width r corr numerically and compare the result with 

r w . 



2. Approach 

Starting point is the expression (see the review [![) 

Sab(E) = S ab - 2inJ2w af ,[D~ 1 (E)]^W ub (3) 

for the element S ao (E) of the scattering matrix connecting channels a and 6, 
with 

{E) = ES^ - +i^Y. w »cW w . (4) 

c 

Here E is the energy. The real and symmetric matrix H with elements 
and n^v = l,...,N is a member of the Gaussian orthogonal ensemble of 
random matrices (GOE). The elements are Gaussian-distributed random 
variables with zero mean values and second moments given by {H^H pa -) = 
(A 2 /iV)[(5 MP (5^ cr + SficrSvp]. The matrix H represents N quasibound levels and 
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their mutual interaction. The parameter A has the dimension energy and defines 
the average level spacing d of the eigenvalues of H . In the center of the GOE 
spectrum we have d = ttX/N. The parameter d defines the energy scale so that 
both E and r corr are expressed in units of d. The real matrix elements W cfJ _ cou- 
ple the space of quasibound levels to A channels labelled a, b, c, In the cases 

considered in the present work the amplitudes for the passage from an intrinsic 
state to a scattering channel coincide with those for the reverse process, that is 
W vc = W cv . Without loss of generality we assume that W all W bfl = Nv 2 5 ab - 
The parameters v\ define the mean strength of the coupling to channel a. Since 
H is random, the S*-matrix is a matrix-valued random process that depends on 
E. All moments and correlation functions of S(E) (defined by averaging over 
the GOE with the energy at or close to the center of the GOE spectrum) depend 
only on the average 5*-matrix elements (S a b), on the transmission coefficients 
T c , and on energy differences. The latter are expressed in units of d. With 
x a = n 2 vl/d we have 

{Sab) = YT^ a Sab ' Ta= or^- (5) 

In Ref. [H , the autocorrelation and cross-correlation functions of the elements 
of the S-matrix are given in terms of these parameters. They are worked out 
for fixed A in the limit N — > oo. We do not repeat the analytical expressions 
here. These contain a threefold integration over real variables. We make use 
of a simplification of these integrals in terms of variable transformations first 



introduced in Ref. 18[ and summarized in the Appendix of Ref. |5j. For a 
given set of transmission coefficients Ti, T2, . . . , Ta the resulting formula for the 
S'-matrix autocorrelation function 

C ab (e) = {S ab (E - e/2)S* ab (E + e/2)) - \(S(E) ab )\ 2 (6) 

is evaluated numerically as a function of e/d. The full width at half maximum 
of that function yields T COIT /d. 



3. Results 

According to the Weisskopf estimate in Eq. ([1} , the correlation width r corr 
should depend only on the number of channels A and on the sum T — ^ c T c 
of the transmission coefficients. To test that assertion, we have for fixed values 
of A and of T with < T < A calculated the width of the autocorrelation 
function Eq. © for several sets of parameters {Ti, T 2 , . . . , 7a}. These are sub- 
ject to the constraints ^ c T c — T and < T c < 1 and were obtained with the 
help of a random-number generator. The number of sets was typically between 
25 and 100. We evaluated C a f,(e), where a and b take either of the channel 
numbers 1 and 2 and found that in the intermediate regime of weakly overlap- 
ping resonances, the widths of Ci2(e), of Cn(e) and of C22(e) vary from set to 
set, in contrast to the above assertion. The deviations of the ratios r corr /r\y 
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from unity are of comparable size in all three cases. Therefore we show in the 
following only results for C\2(e)- 

To test the dependence of T corr on the values of T± , T 2 associated with the 
incident and outgoing channels in the expression for C 12(e) we considered three 
cases. In case I, we chose Ti,T 2 arbitrarily, that is, we did not sort the trans- 
mission coefficients {Ti, T 2 , . . . , T\} by size. In case II (case III) we ordered the 
transmission coefficients such that T\ , T 2 take the maximal values (the minimal 
values, respectively) of all T's. In case II the channels 1 and 2 are the dominant 
ones, in case III they are the most weakly coupled ones. Case II is relevant for 
the microwave experiments described in Refs. [H, M. l20l|. 

In Fig.[T]we consider case I and plot for several values of T given in the panels 
the ratio r corr /r\y versus A for the correlation function C 12(e). The number of 
sets of Ti,T 2 , . . . ,T\ chosen was 25. To test the statistical significance of the 
results we have increased the set size to 100 and did not observe a noticeable 
change. Each set corresponds to a dot in the plot. The dots scatter about 
a mean value that is close to unity. For fixed T (fixed A), the width of the 
cloud of dots decreases with increasing A (increasing T, respectively). The 
width indicates that in contrast to the Weisskopf formula, r corr does depend 
on the values of the individual transmission coefficients. To further test this 
dependence we considered cases II and III. The ratios r corr /rw resulting from 
each of the 25 sets of transmission coefficients form clouds that for both cases 
are very narrow as compared to those shown in Fig. Q] We do not display these 
as they would cover the upper parts of the clouds shown in Fig. [T] in case II, the 
lower parts in case III. 

Figures [2] and [3] serve to quantify these statements. In a plot similar to that 
of Fig.Q] Fig.[2]shows the average of the ratios r corr /rw over the 25 realizations 
versus A for case I (center curve, circles), together with those for case II (upper 
curve, crosses) and case III (lower curve, diamonds). For all values of A and 
T the deviations of the average ratios from unity are largest for case III and 
smallest for case I. In the latter case the average ratio takes values above unity 
for small A and tends to values close to but below unity even for large A. In 
contrast, for case II the average ratio is larger, for case III it is smaller than 
unity for all values of A. In all three cases, the deviations are largest for A values 
between 10 and 20. However, deviations from unity by more than 10 percent 
occur only for T < 8 or so. The curves rapidly tend to unity when A approaches 
T. Then all transmission coefficients take values close to unity. 

Figure [3] shows similarly the root mean square (rms) deviation of the ra- 
tio r corr /r\v from unity (more precisely: the square root of the mean square 
deviation of the ratios from unity). That quantity takes the largest values for 
case III, and is very small unless T < 8. When A approaches T the rms values 
decrease rapidly for all three cases. 

Because of the constraint T < A only few points are at our disposal in the 
regime of largest deviations. For a more thorough test of the Weisskopf estimate 
we, therefore, also considered the case where A < 8 is fixed and T is varied. 
In Fig. [4] we show the averages of the ratios r corr /rw, with all transmission 
coefficients chosen equal, while in Figs. [5] and |6] cases I - III were considered as 
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done above in Fig. [2] Again the deviations of r cor r/rw from unity are largest 
for case III and smallest for case I. The average ratios are smaller than unity 
for all values of T for case III (lower curve, diamonds) while for case I they 
are slightly smaller than unity for small T and eventually reach a value slightly 
above unity when T approaches A, as is also observed in Fig. [5] for comparable 
values of A and T. For case II (upper curve, crosses) the deviations from unity 
are less than 10 percent unless T < 3; the ratio r corr /IV is above unity for 
all T. The dependence of the curves on A is similar to that for equally chosen 
transmission coefficients in Fig. 2] although for given values of A and T the latter 
deviate from unity much less than the former. As in Fig. [5] all curves in Fig. [5] 
tend to unity for T close to A. 

In Fig. [6] we show the rms values for cases I and II. These are less than 0.1 
for all values of T and A that were considered. As suggested by Fig. [5] the rms 
values for case III are always larger than those for cases I and II but remain 
smaller than 0.1 for T > 3. We do not show case III in order not to blur Fig. [SJ 

We also compared the Weisskopf estimate with data from a microwave ex- 
periment described in Refs. 0, [ioj]. The transmission amplitudes S^S^i and 
the reflection amplitudes Sn, S22 of microwaves coupled into and out of a flat 
resonator via two antennas 1,2 were measured. The resonator had the shape 
of a tilted stadium billiard whose classical dynamics is chaotic [2l| . The trans- 
mission coefficients T\ ~ T2 associated with the two antennas were determined 
via the relation T c = 1 — |(5 CC )| 2 from the reflection spectra. Here the angular 
bracket denotes a frequency average. To this end the spectra, which exhibit 
isolated resonances for low frequencies and increasingly overlapping resonances 
with increasing frequency, were devided into 22 frequency intervals of equal size. 
Dissipation into the walls of the resonator was accounted for by introducing 300 
weakly coupled fictitious channels with equal transmission coefficients T c where 
T c <C T\,Ti for c > 3. Hence channels 1 and 2 are the dominant ones in 
these experiments, as in case II. The transmission coeffcients of the absorptive 
channels were determined from a fit of the analytic expression for the autocor- 
relation function C 12(e) to the experimental one 0- The width of the best fit 
yields r corr while Tw is computed from the resulting transmission coefficients 
and Eq. (TTJ . From the measured spectra we thus obtained altogether 22 values 
for the ratios r corr /rw- These are shown in Fig. [7] For T < 1 they are very 
close to unity, as expected for a large number of weakly coupled channels. They 
approach a maximum for T ~ 2 and then decrease again. For T > 4 the de- 
viations of the ratios from unity are around 5 per cent. The ratios are larger 
than unity for all values of T, in agreement with the numerical results for case 
II. Thus, in these experiments T-yv underestimates the correlation width. The 
deviations of the Weisskopf estimate Tw from r corr are similar for the sets of 
transmission coefficients resulting from other microwave experiments described 



in Rcf. 22]. In another experiment [23[ a superconducting chaotic microwave 



billiard was used, and dissipation by the walls was negligible. In that experiment 
three attennas were attached to the resonator so that A = 3. In Ref. [23[ the 
transmission coefficients are determined from the partial widths of the measured 
resonances. Their sum T ~ 0.55 yields Tw while r cor r is determined from the 
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experimental autocorrelation function shown in Fig. 4 of Ref. ■ This yields 
Tcorr/rw — 1.13, in good agreement with our numerical results for small values 
of A and T < A. 

From the results obtained with the non-sorted transmission coefficients we 
conclude that the Weisskopf estimate Eq. ((TJ) constitutes a good approximation 
to r corr for practically all values of ^ C T C . Maximal deviations occur for small 
values of A and T unless A ~ T. That statement applies also when incident 
and outgoing channel are the dominant ones (case II). The largest deviations 
are observed when these are the most weakly coupled channels (case III). In 
cases I and II relative deviations of r corr /rw from unity larger than 10 percent 
are only observed for T < 4. Generally, for T ~ A or for T > 8 the deviations 
are less than 10 percent and decrease rapidly with increasing T or A. We have 
shown that our results are relevant for the microwave experiments described in 
Refs. 0, HO, 22, 23[. It would be interesting to perform similar tests on the data 
obtained, for instance, in the experiments in Refs. 24, 25j |. 
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Figure 1: The ratio r corr /rw (dots) versus the number A of channels for several values of 
T = ^) T c as indicated in the panels. For each value of A and of T, 25 sets of the parameters 
Ti, T2, . . . , 7a were randomly chosen. Each such set corresponds to a dot in the plot. 
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for case II (upper curves, crosses) and for case III (lower curves, diamonds). 
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Figure 7: Ratio r corr /rw for several sets of transmission coefficients taken from a microwave 
experiment 0, l2fj| | versus T. The sets consist of two dominant channels 1 and 2 and 300 
weakly coupled ones, i.e., T c -C T\,T2 for c > 3. 
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